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We measure the transport properties of mechanically strained single crystals of BaFe2(Asi_a;Pi)2 
over a wide range of x. The Neel transition is extremely sensitive to stress and this sensitivity 
increases as optimal doping is approached, even though the transition itself is strongly suppressed. 
Furthermore, we observe significant changes in the superconducting transition temperature with ap- 
plied strain, which mirror changes in the composition x. These experiments are a direct illustration 
of the intimate coupling between different degrees of freedom in iron-based superconductors, re- 
vealing the importance of magneto-elastic coupling to the magnetic and superconducting transition 
temperatures. 
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I. INTRODUCTION 

Materials that exhibit unconventional superconductiv- 
ity are almost always in the proximity of alternative, of- 
ten magnetic ground states. Each ground state is char- 
acterized by a broken symmetry and an associated or- 
der parameter which acquires a finite value at the criti- 
cal temperature, indicating the transition has occurred. 
The iron based superconductors fall within a broad fam- 
ily of correlated electron materials which are related to 
anti-ferromagnetism, joining the cuprate, heavy fermion 
and layered organic superconductors. The relationship 
of the magnetic, structural and other (sometimes un- 
known) order parameters, and particularly how these 
conspire to give rise to high superconducting critical tem- 
peratures Tc, is one of the most important experimental 
challenges in understanding the mechanism behind high- 
temperature superconductivity. 

In the present work, we reveal the intimate relation- 
ship between different broken symmetry ground states in 
the BaFe2As2 superconducting family of iron-pnictides. 
When left chemically unmolested, these materials are 
characterized by high temperature phase that is tetrag- 
onal (Tet) and paramagnetic, transitioning at ~ 138/^ 
to an orthorhombic (Ort), coUinear antiferromagnet 
(AFM)i. In this case, the structural transition breaks 
tetragonal symmetry (C4 C2), and the shear strain 
Uxy = dyUx + dxUy plays the role of the order parame- 
ter. The magnetic order breaks both spin-rotational and 
tetragonal symmetry, characterized by an order param- 
eter corresponding to the staggered sublattice magneti- 
zation Mi where i = 1,2 refers to the magnetization of 
each sublattice^i"— . When BaFe2As2 is electron, hole or 
isovalently 'doped— i^"—, these transitions are suppressed 
and superconductivity (SC) emerges, with optimal ap- 



pearing when magnetism is completely absent, indicating 
that AFM and SC order parameters compete. While for 
electron doped materials the structural (Tg) and mag- 
netic (T/v) transitions separate^iSii^ with Ts > T^, in the 
present isovalently substituted BaFe2(Asi_a;P3;)2 mate- 
rials, no such splitting is observed at any composition. 

Several theoretical descriptions have attempted to ex- 
plain the coupling between the structural and magnetic 
transition, based on pure ferroelastic phenomenology or 
inclusion of a nematic order paramete r^''^'^" . Despite the 
different approaches of each of the works, all of them 
highlight the importance of the magneto-elastic coupling, 
which can cause the two transitions to split or to oc- 
cur simultaneously. Here, we investigate the role of the 
magneto-elastic coupling by studying the thermodynamic 
response of the material BaFe2(Asi_i:Pa;)2 to a small me- 
chanical strain applied along the tetragonal [110]t di- 
rection, or equivalently the orthorhombic b axis when 
T < Tpf. We find that a small shear stress a can signifi- 
cantly alter the Neel transition T/v and superconducting 
Tc, in a manner akin to changing x. Surprisingly, even 
though the magnetism is almost completely suppressed 
at optimal doping, the effect on the magnetic transition 
grows, suggesting that magneto-elastic fluctuations sub- 
stantially increase near optimal Te- 



ll. EXPERIMENTAL METHODS 

The growth of single crystals of BaFe2(Asi_^P2:)2 is 
described elsewhere^i. It is worth noting however that 
the quality of the crystals can be improved by anneal- 
ing within the growth for a week at 900° C. Samples 
were mechanically strained along either the [100] t or 
[110]t direction using a custom built mechanical device 
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FIG. 1. Temperature dependence of the normalised resistivity, /9/p(300if), for unstressed (blue), and uni-axially stressed along 
[110]t (red) BaFe2(Asi_a;P2,)2 crystals with x = 2.6% to a; = 32%. Tn determined from dp/dT are labeled and denoted by 
vertical lines with corresponding colors (see Appendix |X] for details). 



described in RefJ^. A cantilever was pressed against 
the sample by adjusting a turnable screw, applying < 
IQMPa. Each sample was cut to have similar dimen- 
sions (300x200x80/iTO^ - the largest dimension of the 
batch with the smallest samples). To reduce the errors 
associated with differences in pressure, the experiments 
were repeated on 3 samples from each batch. This gives 
us confidence that we are able to apply a similar stress 
for all samples and hence that the changes we detect be- 
tween samples from different batches (different P content 
x) are in fact systematic. This study is distinct from our 
earlier investigations of the transport anisotropy, which 
is a non-equilibrium property, whereas we presently focus 
to the effect of mechanical strain on the temperature of 
the phase transitions themselves. 



III. RESULTS 

Figure[T](a)-(h) illustrates the main experimental data. 
In each panel we show the normalised resistivity vs tem- 
perature for an unstressed (blue) and stressed (red) crys- 
tal at a given doping, where the stress has been ap- 
plied along the [I10]t- The vertical lines denote the as- 
signed Ncel transitions for each curve which have been 
determined by the minimum in the resistivity deriva- 
tive (Appendix Note that in contrast to electron 
doped materials where two anomalies are observed in 
dp/dT in the unstrained samples, we only observe one 
in BaFe2(Asi_a;Pa;)2, indicating that T/v ~ T!s for allas^. 



This may suggest that the magneto-elastic coupling in 
these materials is perhaps largeri^^i^, though we point 
out that other studies have observed split transitions in 
these compounds^. In the presence of mechanical strain, 
the structural transition will naturally broaden to higher 
temperatures, but by continuity, the anomaly seen in the 
data of mechanically strained samples must be associated 
with the magnetic order. 

In all the samples we studied, ISTn = TN{a) — Tisi{a = 
0) > 0, where a indicates the mechanical strain field (the 
stress). Intriguingly, for the unstressed optimally super- 
conducting X = 32% sample (blue curve in Figure [1] (h)) 
there is no detectable magnetic transition, but after ap- 
plication of stress a distinct minimum arises at AhK, al- 
most identical to the minimum seen in the unstressed 
samples at lower doping (consider blue curve in Figure 
[T](g)). It appears that the magnetic transition has been 
summoned from beneath the superconducting dome by 
the application of mechanical strain. Even though the 
magnetic order itself vanishes, the strong magneto-elastic 
coupling as well as magnetic and elastic fluctuations re- 
main. 

In Figure [2] we show the same data presented in Figure 
[U but focus around the superconducting critical temper- 
ature at each composition. In this case Tc is defined as 
the midpoint in the superconducting transition. Even 
though the transition likely gets broader with the appli- 
cation of strain, the superconducting nevertheless can 
be seen to decrease as strain is applied until composi- 
tions beyond optimal x, where we observe Tc to increase 
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FIG. 2. Expanded view of the data shown in Fig [T] around 
to Tc for X = 21.3% to a; = 32% ((a) to (d) .underdoped) and 
X = 21.3 % ((e),overdoped). Red and blue curves correspond 
to the normalized resistivity of unstressed, and uni-axially 
stressed (along [110]t) respectively. For underdoped samples, 
ATc > 0; for overdoped samples, ATc < 0. 
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FIG. 3. Temperature dependence of normalized in-plane re- 
sistivity p/ p(300K), for unstressed crystals and (a) uni-axially 
stressed along [110]t (red), (b) uni-axially stressed along 
[100]t (green), and (c) uni-axially stressed along the crys- 
tal c axis (yellow) BaFe2(Asi-a;Pi:)2 crystals with x=28.7%. 
The upper inserts show dp/dT vs T of each crystal around Tm 
with corresponding colors. For (a), ATiv ~ 21 K, (b), ATat 
~ 5 A", and (c), ATjv ~ —5K. The lower incepts are magni- 
fied plots of /9 vs T around Tc. For (a), ATc 2.5K, (b), 

ATc 0.3K, and (c), ATc ~ ^K. 



(Figure [2](e) illustrates that ATc = Tda) - T^a = 0) 
switches sign beyond optimal x). Finally, the affect on 
T/v and Tc is proportional to the amount of applied 
strain, and this can be demonstrated by applying a sys- 
tematically increasing amount (Appendix [B|. These data 
are suggestive that the application of shear stress has a 
similar affect to decreasing x across the phase diagram. 

To ensure that this effect is indeed intrinsic to pressure 
along one of the orthorhombic axes, we also apply pres- 
sure along [100]t, shown for comparison on two samples 
from the same batch in Figure [3] (a) and (b). Mechanical 
strain in this direction results in small changes in both 



magnetic and superconducting transitions. We further- 
more apply pressure in the inter-layer c direction (Figure 
[3](c)) and found the opposite behavior whereby AT^ < 
and ATc > 0. This suggests that the ratio c/b whether 
directly or indirectly, likely plays a role in the supercon- 
ducting mechanism. 



IV. DISCUSSION 

In Figure |4] (a) we plot the change in the Neel tem- 
perature ATjv as a function of doping for nominally the 
same stress at each doping. Surprisingly, even though the 
magnetic and the structural transitions are suppressed as 
a function of doping, the amount that T/v can be changed 
by stress monotonically increases with P content, within 
our error bars. This is in contrast to the value of the re- 
sistivity anisotropy itself, which is highly non-monotonic 
with doping (Appendix [C]). If we also include data of 
T/v(o') at a; = 32%, ATn appears to have the largest re- 
sponse at optimal doping, as shown in Figure |4] (a). This 
large enhancement of AT/y implies an enhanced suscep- 
tibility of the AFM ground state to shear stress. 

In Figure [4] (b), we illustrate the equivalent plot of 
changes in the superconducting critical temperature AT^, 
as a function of doping for nominally the same strain. 
ATc, unlike AT/v, is not monotonic with x. Compar- 
ing the dependence with the evolution of the unstressed 
superconducting transition with x, it appears that the 
magnitude of ATc is largest where the dome is steepest 
and small otherwise; mathematically ATc oc —dTc/dx. 
Indeed, considering Figure [4] (a) a very similar equation 
likely applies to the magnetic transition, so that AT/v on 
—dTj^/dx. As T/v is increased, there are likely fewer elec- 
trons available to participate in superconductivitj*i^ii^, 
and so the fact that Tc decreases with applied stress in the 
underdoped region is not surprising, since T/v increases. 
However, we note than an unstrained sample with a given 
T/v has always a Tc that is lower than a strained sample 
with the same T/v, i.e. Tc(T/v,0) < Tc(T/v,fT). There- 
fore, there is an intrinsic effect of stress on Tc, beyond 
the indirect effect due to the competition between AFM 
and SC. 

We employ a Ginzburg-Landau (GL) approach to ob- 
tain further insight into our observations. The GL model 
has been applied to the ferro-pnictides by several authors 
already to describe the coupling between structure and 
magnetisrrt^Ti^iiSrJ^. From symmetry considerations, the 
coupling between magnetic and elastic degrees of freedom 
enters the free energy via: 



Fme =9 {Ml- M2) 



(1) 



where g is the magneto-elastic coupling. In order to de- 
scribe the present experiment, we need to add also the 
term —au^y, where a is the applied mechanical stress. 
If one assumes that the structural and magnetic transi- 
tions occur independently, the coupling ([T]) effectively ties 
them together—. Alternatively, it has been proposed that 
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the structural transition is a secondary consequence of 
an underlying electronic order dubbed nematioi^. In this 
case, an independent order parameter 77 oc (Mi • M2) 
condenses and triggers the structural transition via the 
coupling ([1]). Within this approach, the elastic degrees of 
freedom arc not intrinsically soft and can be integrated 
out from the partition function (see Supplemental Ma- 
terial for more details), yielding the contribution to the 
free energy cx ^ (Mi ■ M2), where C° is the bare elas- 
tic shear constant (see also ^^i^° ). This term shows that 
the mechanical stress is converted into a conjugate field 
(Tg/C^ to the electronic order parameter rj. Furthermore, 
it also changes the magnetic part of the free energy, re- 
sulting in an increase of the magnetic transition temper- 
ature TN{a ^0) >Tn{(7 = 0). 

As X approaches optimal compositions, we observe that 
the magneto-clastic response is enhanced, which suggests 
that g/C^ increases substantially and is strongest at opti- 
mal doping. In contrast, our previous mechanical strain 
studies on Ba(Fei_2,Co2^)2As2 did not show significant 
changes in T^r, and changes have only been observed for 
pressures '^5x those used here^, though a recent neu- 
tron study in the parent compound could detect changes 
in Tn at much smaller pressures^^. Nevertheless, as a 
function of Co doping the difference between T/v and T5 
grows^, which could be interpreted in terms of a de- 
creasing g^^. Furthermore, even though the lattice soft- 
ens as a function of temperature, the average value of 
in fact increases as a function of Co doping^i^. Here, in 
contrast, the effect of stress on T/v is strongly enhanced as 
X increases in BaFe2(Asi_a;Pa;)2, suggesting that either g 
becomes larger or (7° smaller, or both. Another possibil- 
ity is that, near optimal doping, where there is no struc- 
tural or magnetic transitions, the nematic susceptibility 
Xnom is enhanced, providing an additional contribution 
that enhances the "conjugate field" f^Xncm (Ml •M2) 
and, consequently, AT^r (see Appendix [D|. Interestingly, 
experiments have indicated that magnetic fluctuations 
are critical at optimally-doped BaFe2(Asi_a;Pa;)2^, which 
could suggest a close connection between nematic and 
magnetic fluctuations in these compounds^ii^. 

Furthermore, magnetic fluctuations enhance the repul- 
sive inter-band pairing interaction that can lead to an 
unconventional superconducting stated. Nematic fluc- 
tuations, on the other hand, give rise to an attractive 
intra-band pairing interaction, which can potentially en- 
hance the transition temperature of the unconventional 
SC phased. 

Previous x-ray studies on Ba(Fei_a;Coa;)2As2 showed 
that Uxy is strongly suppressed below indicating 
that the SC and orthorhombic phases compete. One 
would then expect that by applying stress and inducing 
a non-zero u^y^ Tc would decrease. However, our obser- 
vations that Tc (cr) > Tc (0) in the overdoped region and 
Tc(TN,a) > Tc{Tn,0) in the underdoped region suggest 
that the applied stress may actually lead to an intrinsic 
increase of Tc- To understand the effect of mechanical 
stress on superconductivity, we can compare to the case 
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FIG. 4. (a) r]v(cr = 0) (blue, left axis) and ATjv (pink, 
right axis) vs P concentration x. Points at optimal doping 
are distinguished as open circles because they rely on the 
assumption that riv((T = 0) = 0. (b) Tda = 0) (blue, left 
axis) and ATc (pink, right axis) vs P concentration x. Dotted 
black line indicates ATc = 0. In both (a) and (b) the source of 
error is predominantly related to the our ability to accurately 
determine minima in the resistivity derivative about the Neel 
and superconducting transitions. The size of the effect we see 
on each of these transitions is confirmed on multiple samples. 



of high- Tc copper-oxide based materials^Zr— . Though 
the effects vary between different compounds. Hardy et 
al. proposed a unified picture of the influence of uni- 
axial strain (excluding YBCO), whereby changes in Tc 
could be accounted for by changes in the ratio c/o^S. In 
the present study, c/b will always increase when stress is 
applied along [110]t, and by a smaller amount when ap- 
plied along [100] T, but will decrease when applied along c 
(see Figure [3|) . The changes we are able to invoke on the 
under-doped samples follow this trend, so that qualita- 
tively ATc oc — A(c/6). However, one cannot say whether 
it is the lattice parameters alone, their ratio or some 
other systematically adjusted internal parameter which 
is most important (such as the As-Fe-As bond angle). 
Direct structural measurements as a function of mechan- 
ical strain are required to answer this question. 



V. CONCLUSIONS 

In conclusion, we have found a strongly enhanced 
magneto-elastic response in BaFe2(Asi_a;P2:)2 as x ap- 
proaches optimal doping, which may be related to the 
superconducting mechanism itself. We also find that me- 
chanical strain can directly couple to the superconduct- 
ing order parameter in a manner that is similar to de- 
creasing the P concentration x. These experiments are 
therefore a direct illustration of the subtle coupling be- 
tween different degrees of freedom in BaFe2(Asi_2:Pa;)2. 
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Appendix A: Determination of the Neel 
temperature from the resistivity derivative 

The Fermi surface reconstruction associated with the 
Neel order at Tjv appears as a pronounced minimum in 
the derivative with respect to temperature. After ap- 
plying strain, we observe an increase in T/v, as shown in 
Fig. [5l Blue and red represent unstressed and uni-axially 
stressed (along [ffO]T) I3aFe2(Asi_2;Pa;)2 single crystals 
respectively. Even though the stress broadens the tran- 
sition, the increase in T/v can be easily resolved. 



constants satisfy u > and A > 0, such that in the 
free-energy minimum Mi and M2 are either parallel 
or anti-parallel, corresponding to the two possible mag- 
netic stripe configurations with ordering vectors (tt, 0) 
and (0,7r) in the f-Fe unit call Brillouin zone. It is 
convenient to introduce the order parameters of these 
two magnetic stripe states Ai and A2, such that Ai = 
(Ml + M2) /2 and A2 = (Mi - M2) /2. Notice that 
Mf + M| = 2{Al + Af) and Mi ■ M2 = A? - A^. 

We now consider the nematic part in a phenomeno- 
logical way. The nematic order parameter ip breaks the 
tetragonal symmetry of the lattice. At high enough tem- 
peratures (compared to the structural transition temper- 
ature), such as those for the optimally-doped compounds, 
we consider the free energy expansion only up to second- 
order in the nematic order parameter: 



Appendix B: Systematic response of Tn to pressure 

Although we cannot determine the absolute value of 
stress applied, we can nevertheless tune the amount 
of pressure applied by gradually tightening the screw 
on the device. As a typical example, two samples of 
BaFe2(Asi_:j,P:i,)2, x = 23.1% and x = 28.7% are shown 
in Fig. [HI Each sample was measured with a different 
amount of stress applied along [110]t. IAT/vI and \ATc\ 
increase with increasing stress. 



Appendix C: Resistivity Anisotropy 

The in-plane resistivity anisotropy of 
BaFe2(Asi_a;P3;)2 was obtained by measuring the 
resistivity of mechanically detwined crystals as described 
elsewhere^. Clearly, the resistivity anisotropy has a 
highly non-monotonic dependece on doping. This is in 
contrast to the trend of the response of T/v at constant 
stress as a function of doping, which is a monotonic 
increase with the concentration shown in Figure 21 



Fnen. = ^ (xi°l„) " (A? - A^) (D2) 

where k is the coupling between nematic and magnetic 
degrees of freedom, and Xnem is the static nematic sus- 
ceptibility. For the elastic part, we consider a harmonic 
lattice 



Fel = 'Y'^lv - 9Uxyf ~ Cru^y (D3) 

where g is the magneto-elastic coupling, C° is the bare 
shear modulus, and a is the applied stress. 

To study how T/v^o changes as function of cr, we first 
integrate out the elastic degrees of freedom from the par- 
tition function 



exp 



2C0 



(D4) 



Substituting in Eq. (|D2|1 . the nematic free energy be- 
comes 



Appendix D: Ginzburg-Landau analysis 

To understand how the different degrees of freedom af- 
fect Tn in mechanically stressed samples, we use a phe- 
nomenological Ginzburg-Landau model for the magnetic, 
nematic, and elastic degrees of freedom. For the magnetic 
part, we have: 



Xnem 2 



.{A\ 



A^) 



9o_ 



(D5) 



where we defined the renormalized nematic susceptibility 
Xnem = [xnlm^ - ■ If wc cousidcr the regime 
where the nematic free energy can be approximated by 
the quadratic expansion (|D2p . we can also integrate out 
the nematic degrees of freedom, obtaining 



F„ 



2 



Mt 



Ml) 



[Ml 



Ml) 



X 



(Ml ■ M2)' 
(Dl) 

where Mi and M2 are the staggered magnetization of 
the two interpenetrating Neel sublattices. Here, we de- 
fined tq = a (T — TNfl)i with T/v^o denoting the mean- 
field magnetic transition temperature. The coupling 
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FIG. 5. Derivative of normalised resistivity with respect to temperature, -^p/ p{300K), for unstressed (blue), and uni-axially 
stressed along [110]t (red) BaFe2(Asi_a;Pi)2 crystals with x=2.6% to x=S2%. The minimum in dp/dT is shown in correspond- 
ing color, and is interpreted as the Neel temperature Tm- 




FIG. 6. Temperature dependence of normalized resistivity 
of BaFe2(Asi-^P;r)2, x = 23.1% and x = 28.7%, at three 
systematically increased amounts of stress. The normalized 
resistivity pi with i = {0, crl, a2} stand for unstressed (blue), 
intermediately stressed (green) and highly stressed (red) re- 
spectively. The upper inserts show dpi / dT vs T near Tjv , and 
the lower inserts are expanded plots of pi near Tc. 



FIG. 7. In-plane resistivity anisotropy pal Ph as a function 
of temperature and doping for BaFe2(Asi_a;Pa;)2. The color 
scale has been obtained by a linear interpolation between ad- 
jacent data points. Black circles and squares indicate Ts/Tn 
and Tc respectively. 



F„iag = ra (Al 



-(A? 
2 ^ ^ 



A + 



Substituting in Eq. 
becomes 



Dip , the magnetic free energy 



XnemCg 



(A? - A^) 



{Al-Aiy 



(D7) 
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The last term acts as a conjugate field and breaks the 
tetragonal symmetry, selecting the magnetic stripe con- 
figuration corresponding to the Ai order parameter (or- 
dering vector (tt, 0)). Since = a (T — T/v.o), the mag- 
netic transition temperature is given by 

TN = Tj,,o+^^^^y (D8) 
Hence, the increase in Tn is proportional to the ap- 



plied strain cr. The enhanced response at optimal doping 
can be due to one (or a combination) of the following 
features: an intrinsic softening of the lattice (i.e. de- 
crease C°), an enhancement to the magneto-elastic cou- 
pling (i.e. increase of g and/or ft), and an enhancement 
of nematic fluctuations (i.e. increase of Xnem)- A similar 
enhancement in Tn is also expected even in the absence 
of a nematic order parameter, as pointed out recently by 
Cano and PauU^. 
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